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Aims. The goal of this study is to obtain accurate estimates for the individual masses of the components of the tight binary system T Tau S in 
order to settle the ongoing debate on the nature of TTauSa, a so-called infrared companion. 

Methods. We take advantage of the fact that T Tau S belongs to a triple system composed of two hierarchical orbits to simultaneously analyze 
the motion of T Tau Sb in the rest frames of T Tau Sa and T Tau N. With this method, it is possible to pinpoint the location of the center of mass 
of TTauS and, thereby, to determine individual masses for TTauSa and TTauSb with no prior assumption about the mass/flux ratio of the 
system. This improvement over previous studies of the system results in much better constraints on orbital parameters. 

Results. We find individual masses of 2.73+0.31 M Q for TTauSa and of 0.61±0.17M o for TTauSb (in agreement with its early-M spectral 
type), including the uncertainty on the distance to the system. These are among the most precise estimates of the mass of any Pre-Main 
Sequence star, a remarkable result since this is the first system in which individual masses of T Tauri stars can be determined from astrometry 
only. This model-independent analysis confirms that T Tau Sa is an intermediate-mass star, presumably a very young Herbig Ae star, that may 
possess an almost edge-on disk. 
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1. Introduction 

The most robust method of determining the mass of stars is by 
monitoring the Keplerian orbital motion of multiple systems. 
This method is particularly valuable for contracting Pre-Main 
Sequence (PMS) stars, such as T Tauri stars (TTS), for which 
mass estimates from evolutionary models can differ by factors 
as large as 2 (e.g., Baraffe et al. 2002). However, there are only 
about 20 individual TTS with well-determined masses, all of 
them having radial velocity curves or a Keplerian circumstellar 
disk (see Hillenbrand & White 2004 for a recent review). 

The triple system T Tauri is among the best studied PMS 
multiple system (see, e.g, Beck et al. 2004). The optically 
bright T Tau N, a ~ 2 M Q , ~ 1 Myr-old star, is associated with 
the infrared(IR)-bright T Tau S tight binary system, which was 
resolved less than a decade ago (Koresko 2000). The latter 
consists of TTauSb, a deeply embedded, early M-type TTS, 
and of TTauSa, an extremely red and highly variable object 
(Ghez et al. 1991) that has been classified as an infrared com- 
panion (IRC). Despite two decades of steady observations, the 
nature of that third component remains largely unknown. It 
shares many spectral properties with heavily embedded, Class I 
protostars, but its physical association with two normal TTS 
seems to contradict this hypothesis if the system is coeval 
(Koresko et al. 1997). Alternatively, Ghez et al. (1991) pro- 



posed that TTauSa, which underwent a strong outburst, is a 
low-mass FUOri object. More recently, it has been suggested 
that TTau Sa is an intermediate-mass PMS star surrounded by 
an opaque edge-on circumstellar disk (Duchene et al. 2005). 

Determining the mass of TTauSa from its orbital motion 
within this triple system would help our understanding of its 
exact nature. While it has been suggested that this system may 
be on an unbound orbit (Loinard et al. 2003), other studies have 
found possible closed orbits and derived a total mass of 2-6 M 
for T Tau S (Johnston et al. 2003; Beck et al. 2004). Using a re- 
fined method, Johnston et al. (2004a,b) have further attempted 
for the first time to determine the individual masses of T Tau Sa 
and Sb. All of these studies suggested that the IRC is the most 
massive component of the system. However, they were based 
on simplifying (and sometimes invalid) assumptions regarding 
the flux ratio and/or mass ratio of the T Tau S system, resulting 
in additional sources of uncertainty. 

Here, the combined orbital motion of the triple system is 
revisited using an almost assumption-free, generic approach, 
in order to obtain a precise estimate of the individual stellar 
masses within the TTauS tight system. This Letter is com- 
posed as follows: the fitting method and the observations used 
in this analysis are presented in Section|3 while the results and 
implications are discussed in Section|3] 
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2. Approach and astrometric dataset 

Two types of astrometric data are available on the T Tau mul- 
tiple system. The first type is high-resolution IR data, which 
resolves all three components since 1997. The second type is 
radio data, which has detected for over two decades only two 
of the components, T Tau N and a southern source. Though this 
southern source is shown by Loinard et al. (2003) to be physi- 
cally associated with T Tau Sb, Smith et al. (2003) and Johnston 
et al. (2004a) argue that it might be offset from it by up to C.'Ol- 
0'.'03. In our analysis, we assume that the radio and IR source 
are indeed the same object. Only a much longer IR follow-up 
of the system would allow one to determine whether this as- 
sumption must be dropped. 



Table 1. New astrometric and photometric data for the T Tau 
system. 



UT Date 


Sep. (") 


P.A. O 


AK 




T Tau S 


a-Sb 






1997 Oct. 12 


0.051+0.009 


218+ 


:8 


3.0+0.1 


2004 Dec. 19 


0.116+0.004 


294.1 + 


:1.4 


0.76+0.12 


2005 Nov. 13 


0.119±0.001 


300.6+ 


:1.0 


0.52+0.13 


T Tau N-Sb 


1997 Oct. 12 


0.733+0.009 


181.1+ 


:0.6 


5.6+0.2 


2004 Dec. 19 


0.657±0.001 


193.07+0.17 


3.25+0.06 


2005 Nov. 13 


0.650+0.001 


193.81+0.15 


2.42+0.12 



2.1. Fitting Method 

Two general approaches have been used in previous attempts 
to fit the orbital motion of the T Tau Sa-Sb tight binary system. 
First, it is possible to consider in the fit only the IR datasets, 
which provide the location of both components. In this case, 
however, the short timespan probed since 1997 poorly con- 
strains the orbit and a wide range of total system masses is 
allowed (Beck et al. 2004). Including our most recent IR mea- 
surements (see Section l2~2l only slightly improves this conclu- 
sion. 

In a second approach, the separation between T Tau N and 
T Tau Sb, estimated from the much longer radio observations, 
is transformed into a TTau Sa-Sb separation prior to fitting the 
orbital motion. In this case, one has to make some assump- 
tions regarding the unknown location of TTau Sa with respect 
to TTauN. In the past, two main categories of assumptions 
were made. First, some authors considered that TTauSa fol- 
lows a constant linear motion with respect to T Tau N (Loinard 
et al. 2003; Johnston et al. 2003), implicitly assuming that it is 
coincident with the center of mass of TTau S. However, while 
TTau Sa is indeed more massive than its companion, the mass 
of the latter is not negligible (see Section^. Second, it was as- 
sumed that the historical location of the photocenter of T Tau S, 
which has been observed for over 15 years, can be used to trace 
the location of T Tau Sa prior to 1997 (Johnston et al. 2004a,b). 
This requires one to assume that the IR flux ratio within the 
T Tau S system remains constant. Evidence for large variability 
in at least one of the two components of the system is growing, 
however (Ghez et al. 1991; Beck et al. 2004). 

To limit as much as possible the number of assumptions 
in the analysis, we decided to fit simultaneously two types of 
data, namely the location of TTau Sb in two rest frames: with 
respect to TTauSa and to TTauN. We did not use the sepa- 
ration between TTauN and TTauSa because it would not be 
independent from the other two types of measurements. 

We assume that TTauSb follows a bound Keplerian or- 
bit around TTauSa and use the 141. 5+2. 8 pc distance to the 
system estimated by Loinard et al. (2005). We further assume 
that the motion of the center of mass of TTauS, which is 
a Keplerian orbit around T Tau N, can be approximated by a 
parabolic motion over the 20 years of observations probed here. 
This is justified as the period of this outer orbit is at least 400 



years (Johnston et al. 2004a,b). Fitting the motion of TTauSb 
in both rest frames simultaneously allows us to determine in- 
dividual masses from astrometry only, as the ratio of the am- 
plitude of the motion of T Tau Sb in both frames is a direct 
measurement of the mass ratio. 

We use a Levenberg-Marquardt x 2 minimization routine 
(Press et al. 1992) that simultaneously fits for 14 independent 
parameters. The TTau Sa-Sb orbit is described by 7 parame- 
ters: semi-major axis (a), eccentricity (e), inclination with re- 
spect to the plane of the sky (/), position angle of ascending 
node (Q, measured Eastward from North), azimuth of peri- 
astron (oj), time of periastron passage (tp) and orbital period 
(P). The parabolic motion of the center of mass of TTauS 
with respect to T Tau N is described by 6 parameters defined 
by its right ascension = oq + v a (t - to) + a a (t - to) 2 and a 
corresponding equation for its declination, with to = 2000.0. 
Finally, the last parameter describes the mass ratio in the sys- 
tem: ii = Msb/Ms, where M sb and M s are the mass of T Tau Sb 
and of the whole TTauS system, respectively. Kepler's third 
law is used to determine Ms from the values of a and P; the 
individual stellar masses Ms a and Msb are derived from the val- 
ues of Ms and /j. The minimization routine was run with several 
tens of thousand initial guesses spanning wide ranges for all fit 
parameters to ensure convergence on the absolute minimum of 
the x 2 function. 

The uncertainties on the fit parameters are estimated us- 
ing a bootstrap method, a robust method when performing a fit 
with many (correlated) parameters. We created 10000 artificial 
datasets containing as many points as the observed dataset, in 
which each point is randomly drawn from a Gaussian distri- 
bution centered on the actual measurement and whose width 
is given by the associated astrometric uncertainty, and run the 
X 2 minimization routine for each realization. The resulting dis- 
tribution of values for each fit parameter is characterized by a 
mean value equal to that of the best fit and by a standard devi- 
ation that we consider as the estimate of the uncertainty on the 
fit parameter. 

Because this fitting method is based on a fit of all available 
measurements without further transformation, and because no 
parameter is set to an ad hoc fixed value, it is very robust and 
essentially assumption-free. 
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2.2. IR astrometric dataset 

A dozen high spatial resolution IR observations of the TTau 
triple system have been published; they usually provide simul- 
taneously the position of T Tau Sb with respect to both other 
stars. All astrometric measurements listed in Beck et al. (2004) 
and Duchene et al. (2005) are included in our fit 1 . Two addi- 
tional observations of the system were recently published by 
Mayama et al. (2006). However, they are significantly offset 
(> 3cf) from all other quasi-contemporaneous measurements; 
these measurements are likely biased or suffer from underes- 
timated uncertainties and they have been discarded from the 
analysis. 

In addition to these measurements, we present new TT-band 
speckle interferometric datasets obtained on the 10m Keck I 
telescope in 1997, 2004 and 2005. For the last two epochs, the 
TTauSa-Sb system is sufficiently well-resolved that the rel- 
ative astrometry of all components could be determined from 
shift-and-add images. On the other hand, in 1997, the system 
was very tight (» 1 .25 A/ D) and traditional speckle analysis was 
conducted. Usual basic reduction techniques were first applied 
to all datasets (see Ghez et al. 1993 for more details). Because 
T Tau is a hierarchical triple system, a 2-step method identical 
to that described in Duchene et al. (2003) was used to extract 
its properties. The new astrometric and photometric data are 
presented in Tabled 

2.3. Radio astrometric dataset 

Loinard et al. (2003) and Johnston et al. (2003, 2004b) have 
published a series of VLA centimeter observations of TTau 
spanning 20 years, from 1983 to 2003. Both groups indepen- 
dently analyzed the 1983, 1988, 1992, 1995 and 2001 datasets, 
and obtained very similar results. We found that using either 
sets of results for these epochs yields very similar fits 2 . In the 
following, we use the 2 cm astrometric results from Johnston 
et al. (2003, 2004b) in order to ensure a maximum uniformity 
among the radio datasets and their analysis, as well as the 1998 
3.7 cm measurement from Loinard et al. (2003), since no 2 cm 
observations were conducted at that epoch. 



3. Results and implications 

Overall, we considered in our fit a total of 1 1 astrometric mea- 
surements of the TTauSa-Sb system and 16 measurements of 
the TTauN-Sb system. This is the most extensive dataset fit- 
ted to date on this system, considering that we did not try to 
fit the historical location of the photocenter of T Tau S (but see 
SectionlTTll. 



1 The astrometric measurements presented in Duchene et al. (2002) 
suffered from an incorrect astrometric calibration: the position angles 
used in this study are 1.1° larger than was previously reported. 

2 It is however incorrect to simultaneously fit both, as was done by 
Tamazian (2004), since they come from the same raw datasets. 



Table 2. Parameters of the best orbital solution. 



Fit parameters 



a(10- 3 ") 


82.1 


+ 1.8 


<*o (10- j ") 


-20.7 


±4.0 


e 


0.466 


±0.031 


£o(10- 3 ") 


-694.3 


±2.9 


in 


37.2 


±5.1 


v ff (10- 3 7yr) 


-9.12 


±0.54 


no 


115.2 


±6.8 


v s (10- 3 "/yr) 


1.71 


±0.39 


w(°) 


12.5 


±10.3 


a a (10- 4 "/yr 2 ) 
a 5 (10- 4 7yr 2 ) 


1.09 


±0.50 


tp 


1995.86 


±0.21 


1.18 


±0.58 


P(yr) 


21.66 


±0.93 




0.183 


± 0.044 


Derived quantities 2 


A (AU) 


11.59 


±0.35 


M s (M Q ) 


3.34 


±0.36 


M Sa (Mq) 


2.73 


±0.31 


M sb (M ) 


0.61 


±0.17 



a The uncertainty on the derived quantities include that induced by 
the distance estimate. 



3.1. Individual masses within TTauS 

Fitting all measurements described in Section|2j a best fit so- 
lution with;^ 2 ed = 1.88 is obtained, a slightly large value con- 
sidering that there are 40 degrees of freedom. In fact, the 1998 
VLA measurement reported by Loinard et al. (2003) is 4.1cr 
away from the fitted value and contributes a large amount to 
the total x 1 - All other measurements are within 3<x of the fit, 
and all but three less than 2<x away from it, as expected from 
a Gaussian distribution of the astrometric uncertainties. This 
particular VLA measurement is the only one included in our fit 
that was obtained at 3.7cm. Substantial differences were noted 
in the past between 2cm and 3.7cm VLA measurements for this 
source (Johnston et al. 2003, 2004b), and we therefore consider 
this point as suspicious. 

Excluding the 1998 VLA measurement, a best fit solution 
which is less than 2.7cr away from all remaining measurements 
and that has a much more acceptable — 1-37 is obtained 
with 38 degrees of freedom. No fit parameter is different by 
more than 1.4<x from the previous solution, and we consider 
this new solution as our best fit in the following. The best fit 
parameters are listed in Tableland the solution is illustrated in 
Fig-H Note that it is not possible to distinguish the ascending 
and descending nodes of the orbit from astrometry only: an 
identical projected orbit could be described with Q = 295.2° 
and to = -167.5°. 

We then performed a consistency check regarding the his- 
torical measurements of the IR photocenter of the T Tau S bi- 
nary. While strong variability prevents using these measure- 
ments in the fit, it must be verified that the photocenter is lo- 
cated along the line joining the fitted location of TTau Sa and 
T Tau Sb at all epochs. Satisfyingly, the photocenter lies within 
2.4<t of the predicted line connecting the two components for 
all measurements listed in Ghez et al. (1995) and Roddier et al. 
(2000), giving more confidence in the validity of the fit. 

Including the uncertainty induced by the distance estimate, 
the total relative uncertainty on Ms a is about 11%, among the 
most accurate individual mass estimates for a TTS that does 
not belong to an eclipsing binary system (Hillenbrand & White 
2004). The uncertainty on Msb is about 28%, still a remarkably 
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Fig. 1. Observed motion of T Tau Sb in the reference frame of T Tau Sa (left) and T Tau N (right). The dotted curves represent 
our best combined fit; each observation is connected to its modeled counterpart by a straight line. In the left-hand panel, the 
square indicates the fixed location of TTauSa and the dashed line connects it to the periastron. In the right-hand panel, filled 
diamonds represent the VLA measurements; the open diamonds indicate the 1998 VLA measurement excluded from the final fit 
and its fitted counterpart. In this panel, the dashed curve represents the motion of T Tau Sa. 



accurate estimate for a TTS, considering that only astrometric 
measurements were taken into account in this analysis. 



3.2. Discussion 

The results obtained here are in good agreement with those pre- 
sented by Johnston et al. (2004a,b), though with uncertainties 
that are at least twice as small for all parameters: the best esti- 
mate for Msa to date was 2. 1 +0.8 M e , for instance (Johnston et 
al. 2004b). This significant improvement is due to a combina- 
tion of factors: the increased time coverage of high-resolution 
IR datasets, a simpler characterization of the T Tau N-S orbit 
(parabolic motion instead of a Keplerian orbit with several pa- 
rameters fixed to ad hoc values), the improved distance esti- 
mate obtained by Loinard et al. (2005), and the smaller num- 
ber of assumptions used in our analysis. In this work, the only 
assumptions that could be disputed are the parabolic approx- 
imation of the orbital motion of T Tau S around T Tau N over 
20 years and the match between the radio and IR location of 
T Tau Sb. The fact that a very satisfying fit was obtained ap- 
pears to validate these hypotheses. 

The mass derived for TTauSb is in excellent agreement 
with its known, early-M spectral type (Duchene et al. 2002), 
as already discussed by Johnston et al. (2004a,b). For the first 
time, however, an accurate estimate of the mass of TTauSa 
is obtained, leaving little doubt that it is somewhat more mas- 
sive than T Tau N itself and suggesting that T Tau Sa is a young 
embedded HerbigAe star. This is the first time that this fact 
is established with such certainty for an IRC. As discussed in 
Duchene et al. (2005), it is likely that TTauSa is surrounded 
by a compact and partially opaque circumstellar disk that ac- 
counts for the apparent faintness of this intrinsically bright ob- 
ject. Given the orbital parameters of the TTauSa-Sb system, 
the radius of this disk must be < 3 AU, i.e., half the periastron 
distance. We remark that this disk must be very inclined with 



respect to the orbital plane (A/ > 45°), raising the possibility of 
gravitational instabilities. 

Within 10 years or so, spatially-resolved IR observations 
of the T Tau triple system will allow the derivation of accurate 
orbital parameters for T Tau S without making use of the radio 
datasets. With such an analysis, it will be possible to deter- 
mine whether the southern radio component of T Tau is indeed 
arising from the chromosphere of T Tau Sb or whether it is a 
separate object that is simply associated with this object. 

On the shorter term, the dynamical stability of the TTau 
system must be studied in more detail: the system contains 
three stars within 100 AU of each other, each of them possess- 
ing a circumstellar disk. The overall velocity of T Tau S with 
respect to T Tau N is on order 6 km/s, very similar to the ve- 
locity of a 100 AU -radius circular orbit for a 5 M© system (as- 
suming Mn = 2M Q ). This suggests that the outer orbit is seen 
close to face-on and has a low eccentricity. If this is indeed 
the case, the application of the stability criterion proposed by 
Eggleton & Kiseleva (1995) suggests that the system could be 
marginally stable. Moreover, the stability of the circumstellar 
disks is not necessarily ensured in such a tight system. This mo- 
tivates a follow-up dynamical study similar to that performed 
for the quadruple system GG Tau (Beust & Dutrey 2005, 2006) 
to determine the possible solutions for the outer orbit that pre- 
serve the three components and their disks for at least 1 Myr, 
i.e., the age of the system. A preliminary analysis, which will 
be presented in a forthcoming paper, suggests that the three- 
body system is stable over a few Myrs, whereas the fate of the 
circumstellar disks is not yet established. 
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